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ABSTRACT 

Aims. We present Herschel/HIFI observations of fourteen water lines in W43-MM1, a massive protostellar object in the luminous star 
cluster-forming region W43. We aim at placing this study in the more general context of high-mass star formation. The dynamics of 
such regions may either be the monolithic collapse of a turbulent core, or competitive accretion. Water turns out to be a particularly 
good tracer of the structure and kinematics of the inner regions, allowing an improved description of the physical structure of the 
massive protostar W43-MM1 and an estimation of the amount of water around it. 

Methods. We analyze the gas dynamics from the line profiles using Herschel-HIFI observations (WISH-KP) of fourteen far-IR water 
lines (Hj^O, H'^O, H'^^O), CS(ll-lO), and C'**0(9-8) lines, and using our modeling of the continuum spectral energy distribution. 
The spectral modeling tools allow us to estimate outflow, infall and turbulent velocities, and molecular abundances. We compare our 
results to previous studies of low-, intermediate-, and other high-mass objects. 

Results. As for lower mass protostellar objects, the molecular line profiles are a mix of emission and absorption, and can be de- 
composed into 'medium' (FWHM^5-10 km s"'), and 'broad' velocity components (FWHM=: 20-35 km s"'). The broad component is 
the outflow associated with protostars of all masses. Our modeling shows that the remainder of the water profiles can be well fitted 
by an infalling and passively heated envelope, with highly supersonic turbulence varying from 2.2 km s"' in the inner region to 3.5 
kms"' in the outer envelope. Also, W43-MM1 has a high accretion rate, between 4.0 x 10"* and 4.0 x 10"^ Mgyr"', derived from 
the fast (0.4-2.9 km s" ' ) infall observed. We estimate a lower mass limit of gaseous water of 0. 1 1 Mq and total water luminosity of 1 .5 
Lq (in the 14 lines presented here). The central hot core is detected with a water abundance of 1.4 xlO"'' while the water abundance 
for the outer envelope is 8 xlO"^. The latter value is higher than in other sources, most likely related to the high turbulence and the 
micro-shocks created by its dissipation. 

Conclusions. Examining water lines of various energies, we find that the turbulent velocity increases with the distance to the cen- 
ter While not in clear disagreement with the competitive accretion scenario, this behavior is predicted by the turbulent core model. 
Moreover, the estimated accretion rate is high enough to overcome the expected radiation pressure. 

Key words. ISM: molecules - ISM: abundances ~ Stars: formation - Stars: protostars - Stars: early-type ~ Line: water profiles 



1. Introduction 

After Hydrogen and Helium, Oxygen is the most abundant el- 
ement in the Universe and water is one of the most abundant 
molecules. Due to the water in the Earth's atmosphere, space- 
based observations are necessary to study the vast majority of 
the water transitions. Because of its importance for life and its 
sensitivity to dynamical, thermal, and chemical processes in the 
interstellar medium, the physics and chemistry of water is one 
of the main drivers of the Herschel Space Observatory mission 
(hereafter Herschel, Pilbratt et al.||2010[l and particula rly of the 
HIFI spectroscopy instrument ( |de Graauw et al.|2010) l. 

Massive stars are rare but are the major contributors to the 
matter cycle in the Universe due to their short lifetimes and rapid 
ejection of enriched material. OB stars dominate the energy bud- 



* Herschel is an ES A space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- 
tant participation from NASA. 



get of star-forming galaxies and are visible at great distances. 
Their formation, however, is not well understood and the clas- 
sical scheme for low-mass star formation cannot be applied as 
such to OB stars. Indeed, young OB stars and protostars strongly 
interact with the surrounding massive clouds and cores, leading 
to a complex and still not clearly defined evolutionary sequence. 
'We generally identify in this sequence objects ranging from mas- 
sive pre-stellar cores, high-mass protostellar objects (HMPOs), 
hot molecular cores (HMC), and finally to the more evolved 
Ultra Compact HII regions stage where the central object be- 
gins to ionize the surrounding gas (e.g., |Beuther et al. 2007 ). The 
classification adopted above is not unique but is consistent with 
the analysis of massive young stellar objects and HII regions in 
the Galaxy made by 'Mottram et al. ( 201 l| l. The most problem- 
atic issue in understanding the massive star formation process is 
to explain how to accumulate a large amount of mass infalling 
within a single entity despite radiation pressure. Models consid- 
ering a protostar-disk system (e.g., Yorke & Sonnhalter 2002 
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|Krumholz et al.|2005[ |Banerjee & Pudritz|2007| l now quite suc- 
cessfully address how the accretion of matter overcomes radia- 
tive pressure. Actually, two main theoretical scenarios have been 
proposed to form high-mass stars, both requiring the presence of 
a disk and high accretion rates: (a) a turbulent core model with 
a monolithic collapse scenario ( |Tan & McKee|2002 McKee & 



estimates toward four other massive YSOs have been made by 



Marseille et al. (2010b I using the lowest two P-H2O lines com- 



|Tan| l2003 ); (b) a highly dynamical competitive accretion model 
involving the formation of a cluster (Bo nnell & Bate 2006). In 
ionized HII regions, star formation could also be triggered by 
successive generations of stars (e.g., |Deharveng et al. 20091). One 
of the implications of the turbulent core model is that molecular 
line profiles should show supersonic turbulent widths while the 
competitive accretion produces cores which are subsonic, hence 
with rather small turbulence (see |Krumholz & Bonnell|2009| l. In 
the deeply embedded phase of star formation, it is only possi- 
ble to trace the dynamics of gas through resolved emission-line 
profiles, such as obtained with HIFI. 

The Guaranteed-Time Key Program Water In Star-forming 
regions with Herschel (WISH, van Dishoeck et al.|201 \) probes 
massive star formation through water observations using the 
HIFI and PACS (P oglitsch et al.|,20 10) instruments. In particu- 
lar, the dynamics of the central regions are characterized through 
the water lines and the amount of cooling they provide is mea- 
sured. In molecular clouds, water is mostly found as ice on dust 
grains, but at temperatures T >100 K the ice evaporates, increas- 
ing the gas-phase water abundance by several orders of mag- 
nitude ( |Fraser et aL]|2001| [Aikawa et aL]|2008| l. This suggests 
that water emission almost exclusively probes the warm inner 
regions of protostellar objects. In order to collapse, the gas must 
be able, among several other conditions, to release enough ther- 
mal energy; a major WISH goal is to determine how much of the 
cooling of the warm region {T >100 K) is due to H2O. 

The high velocity resolution of the HIFI instrument allows us 
to study the dynamics of the gas, detect outflows, and estimate 
infall and turbulent velocities present in protostellar envelopes. 
Even though the water abundance in the envelope is small, low 
energy lines are highly absorbed by the envelope, producing a 
mixture of emission and absorption that requires the spectral res- 
olution of HIFI in order to be studied. The high-energy lines are 
observed in emission and directly probe the warm regions since 
their lower energy levels are not excited in the envelope, and as 
such cannot absorb emission from the hot core. Thus, observa- 
tions of a mixture of low- and high- lying lines of water and its 
isotopologues allow a 'tomography' of the structure of protostel- 
lar envelopes. 

A first review of the WISH KP and results have been pre- 
sented and compiled by [van Dishoeck et aL ( 2011[ l for a few 
objects only. We expect to confirm these results with the study 
of the whole sample. Interestingly, the water line profiles (tran- 
sitions P-H2O In-Ooo and P-H2O 2o2-lii) obtained in a low- 
(NGC1333,|K ristensen et al.|2010| l, an intermediate- ( NGC7129, 



Johnstone et al. 2010), and a high-mass (W3IRS5, [Chavarriaj 



et al. 2010 1 young stellar object (hereafter YSO) exhibit sim- 
ilar velocity components: a broad (FWHM ~ 25 kms ') and a 
medium (~5-10 km s"') one. A naiTower component (<5 km s"') 
is also observed, but not in the massive object, although it has 
been detected toward W3IRS5 by Wampfler et al. (201 1 1 in the 



OH emission. On the opposite, the water emissions differ in in- 
tensity, the massive objects being the strongest emitters. Another 
high-mass star forming region, DR21(Main) has been studied by 



van der Tak et al.| ( |2010j ) in the para ground-state water and '^'CO 
lines. They derived a very low water abundance of a few 10"^ 
in the outer envelope while it is enhanced by nearly three or- 
ders of magnitude in the outflow (7 xlO"^). Water abundance 



bined with the ground state p-Hj^O line. Variations of the outer 
envelope water abundance are measured (from 5 xlO"'" to 4 
xlO but are not correlated with the luminosity or evolution- 
ary stage. Finally, no infall has been firmly established so far for 
these objects. 

In this paper we analyze the water observations towards 
the massive dense core W43-MM1 in the mini-starburst region 
W43, located near the end of the Galactic bar at a distance of 5.5 
kpc ( |Motte et alTp()03 ; Nguyen Luong et al.|201 \\ . This region 
is one of the most interesting star formation regions in several 
aspects (e.g., |Bally et aL]|2010| l. It has a star-forming rate that 
has increased by an order of magnitude compared to ~ 10* years 
ago (Nguyen Luong et al. 201 l| l. The source is an IR-quiet dense 
core following the definition of |Motte et al. ( |2007| l with bolo- 
metric luminosity 2.3 xlO"^ Lq and size of 0.25 pc. This mas- 
sive dense core is expected to be fragmented on small scales and 
should be considered as a protocluster. W43-MM1 is the most 
massive dense core (M~ 3600 Mq, Mott e et aLl|2003| ) of the 
four identified main millimeter fragments of which only MMl is 
within the HIFI beam. The mid-IR emission from the hot core is 
undetected because it is totally absorbed by the envelope, while 
the source is detected at 24 yum (peak flux= 0.3 Jy). A methanol 
and a water maser are detected but there is no near-infrared or 
centimeter emission. Infall (up to 2.9 kms"') has been identi- 
fied by [Herpin efaT] p009 ) from CS data. W43-MM1 resem- 
bles a low-mass Class protostar, i.e., a YSO in its main ac- 
cretion phase, but scaled up in mass and luminosity. However, 
although it looks like a very young mid-IR quiet HMPO, W43- 
MMl has already developed a hot core with temperatures higher 
than 200 K ( [Motte et al.|2003[ [Marseille et al.|2010b[ ), making 
this object very peculiar. All these characteristics make W43- 
MMl a promising massive object to study, potentially rich in 
water: among first 10 massive sources, from different evolution- 
ary stages, observed by [Marseille et al.[ ( [2010a| l, W43-MM1 is 
one of the three only sources (with DR210H and IRAS18089- 
1732) to be detected in the p-H^^O 3i3-32o line. 

We present here new Herschel-HIFI observations of fourteen 
far-IR water lines (Hj^O, Hj^O, Hj^O), providing the most com- 
plete coverage of the water energyladder to date, plus CS(11- 
10) and C'^0(9-8) data. Sections |2] and [3] present our observa- 
tions and results, whose analysis is given in Sect. |4] Then we 
model the observations using a radiative transfer code in Sect. 
[5] We estimate outflow and infall velocities, turbulent velocity, 
molecular abundances, and the physical structure of the source. 
Finally, we discuss (Sect. |6]l the results in terms of massive star 
formation scenarios. 



2. Observations 

Fourteen water Unes as well as the CS(ll-lO) and C'^0(9-8) 
lines (see Table 1) have been observed with HIFI at frequencies 
between 538 and 1670 GHz towards W43-MM1 in March, April, 
and October, 2010 (OD 293, 295, 310, 312, 333, 338, 339, and 
531). The position observed corresponds to the peak of the mm 
continuum emission from Motte et al.[ ( |2003l ) (RA= 18:47:47.0, 
DEC=-01:54:28 J2000). The observations are part of the WISH 
GT-KP 

Data were taken simultaneously in H and V polarizations 
using both the acousto-optical Wide-Band Spectrometer (WBS) 
with 1 . 1 MHz resolution and the digital auto-correlator or High- 
Resolution Spectrometer (HRS) providing higher spectral reso- 
lution. We used the Double Beam Switch observing mode with a 
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Table 1. Herschel/HIFI observed line transitions in W43-MM 1 . Frequencies are from Pearson et al. ( 1 99 1 
Roelfsema et al.| ( [2012 1. The rms is the noise at the given spectral resolution 6v. The energy ot the upper 

- -- - -'"o. 



be the same for Hj'O and 



Beam and t]^^ are from 
level. Ell, is considered to 



Water species 


Frequency 


Wavelength 


E, 


HIFI 


Beam 


''mb 




6v 




rms 


obsid 




[GHz] 


Ipm] 


[K] 


band 


["] 




[K] 


[MHz] 


[s] 


[mK] 




o-Hf0 1,o-loi" 


547.6764 


547.4 


60.5 


la 


38.7 


0.754 


80 


0.24 


1462 


43 


1342219193 


o-Hf O l,o-loi 


552.0209 


543.1 


61.0 


la 


38.0 


0.754 


71 


0.24 


158 


46 


1342192345 


p-Hf0 2o2-lii 


994.6751 


301.4 


100.6 


4a 


21.3 


0.741 


288 


0.12 


502 


107 


1342194990 


o-Hf 3i2-3o3 


1095.6274 


273.8 


248.7 


4b 


19.2 


0.737 


379 


0.24 


1723 


54 


1342194806 


p-HfO 1,1 -Ooo 


1101.6982 


272.1 


52.9 


4b 


19.2 


0.736 


393 


1.1 


3566 


17 


1342191670,1342207372 


p-Hf l„-Ooo 


1107.1669 


272.1 


52.9 


4b 


19.1 


0.737 


379 


0.24 


1723 


29 


1342194806 


o-H^'0 2,2-loi 


1662.4644 


180.3 


113.6 


6b 


12.8 


0.708 


1407 


1.1 


1212 


176 


1342192575 


O-H2O lio-loi" 


556.9361 


538.3 


61.0 


la 


38.0 


0.754 


80 


0.24 


1462 


43 


1342219193 


P-H2O 2ii-2o2 


752.0332 


398.6 


136.9 


2b 


28.2 


0.749 


88 


0.12 


262 


88 


1342194565 


P-H2O 2o2-lll 


987.9268 


303.5 


100.8 


4a 


21.3 


0.741 


340 


1.1 


418 


67 


1342191616 


O-H2O 3l2-3o3 


1097.3651 


273.2 


249.4 


4b 


19.2 


0.737 


379 


0.24 


1723 


54 


1342194806 


P-H2O ill -Ooo 


1113.3430 


269.0 


53.4 


4b 


19.0 


0.736 


393 


1.1 


3566 


17 


1342191670,1342207372 


0-H2O 221 -2i2 


1661.0076 


180.5 


194.1 


6b 


12.8 


0.708 


1407 


1.1 


1212 


225 


1342192575 


O-H2O 2i2-loi 


1669.9048 


179.5 


114.4 


6b 


12.7 


0.708 


1407 


1.1 


1212 


249 


1342192575 


CS(11- 10) 


538.6888 


556.5 


155.1 


la 


38.0 


0.754 


71 


0.24 


158 


46 


1342192345 


C"*0(9-8) 


987.5604 


303.6 


237.0 


4a 


21.3 


0.741 


340 


0.5 


418 


54 


1342191616 



Notes. " This line was mapped in OTF mode. 



throw of 3'. HIFI receivers are double sideband with a sideband 
ratio close to unity. The frequencies, energy of the upper lev- 
els, system temperatures, integration times and rms noise level 
at a given spectral resolution for each of the lines are provided 
in Table 1. Calibration of the raw data into the Ta scale was 
performed by the in-orbit system p ioelfsema et al.|2012 |; con- 
version to Tmb was done using a beam efficiency given in Table[T] 
and a forward efficiency of 0.96. The flux scale accuracy is es- 
timated to be between 10% for bands 1 and 2, 15% for band 3, 
4, and 20 % in bands 6 and 7. Data calibration was performed 
in the Herschel Interactive Processing Environment (HIPE, Ott] 
2010 1 version 6.0. Further analysis was done within the CLASSTj 



package. These lines are not expected to be polarized, thus, after 
inspection, data from the two polarizations were averaged to- 
gether For all observations, eventual contamination from lines 
in the image sideband of the receiver has been checked and none 
was found. Because HIFI is operating in double-sideband, the 
measured continuum level (in the Figures and the Tables) has 
been divided by a factor of 2. 



3. Results 

The spectra including continuum emission are shown in Figures 
[1] and |2] for the rare isotopologues (H^O, Hj**©) and for 
H^*'0 respectively. In addition. Fig |3] displays the CS (11 - 10) 
and C'^O (9 - 8) spectra. We show the HRS spectra, except for 
the p-Hi'^O In-Ooo, o-Hi^O 2i2-loi, P-H2O 2o2-lii, p-HjO In- 
Ooo, 0-H2O 22i-2i2, and 0-H2O 2i2-loihnes where WBS spectra 
were used since the velocity range covered by the HRS was not 
sufficient. 

We derive the peak main beam temperatures and half 
power line-widths for the different line components from multi- 
component Gaussian fits, made with the CLASS software (line 
parameters are given in Table |2]i. All lines associated with the 
source have V/.„. ^ 95 - 100 kms"'. Since W43-MM1 is located 
in a very crowded area of the galactic plane (b~ -0.1°), fore- 



ground clouds contribute to the spectra through water absorp- 
tions at Wisr shifted with respect to W43 velocity in the 0-H2O 
lio-loi> P-H2O lii-Ooo and 0-H2O 2i2-loi fines spectra. These 
will be discussed in Appendi x [A] 

As in [Johnstone et aL] ( |2010| ), [Kristensen et al.| ( |2010| l. 



and Chavarria et al. (20101, we adopt the following terminol- 
ogy for the different velocity components: broad (FWH]VI^20- 
35 kms"'), medium (FWH]VI^5-10 kms"^), and narrow (~3 
kms"'). 



3.1. C-bearing species 

Although the main goal of WISH is to observe H2O, the HIFI 
bands included the C'^O (9-8) and CS (11-10) fines as weU. 
They are seen in emission without any broad component contri- 
bution. While the source velocity derived by IVIotte et al. ( 2003| l 
from H'^CO^ observations is 98.8 kms the best Gaussian fit 
for the observed C'**0 (9 - 8) line emission is centered at 99.4 
kms"' (with a FWHM of 5.6 kms"'), which is consistent tak- 
ing into account the uncertainties. The CS 11-10 line profile is 
best fit by two Gaussians of 2.9 and 10.2 kms"'width, centered 
at 98.4 and 98.8 kms"', respectively, but the CS line profile is 
asymmetric, possibly revealing infall, blue-shifting the emission 
peak. In the following, we adopt V/s^ - 99.4 km s"'as the source 
velocity. 

Note that an unidentified line is detected in absorption in the 



http://www.iram.fr/IRAMFR/GILDAS/ 



C^^O (9-8) spectra at 88.6 kms"', i.e., 987.594/972.747 GHz in 
USB/LSB). 



3.2. HfO lines 

Most of the Hj^O lines can be described as the sum of a medium 
(FWHM^5-10 kms"') and a broad (up to 35 kms"') velocity 
components. The broad one is likely the signature of an out- 
flow. The observed dependence of the outflow velocity extent 
with Eyp is probably due to sensitivity effects (the outflow being 
more difficult to detect with weaker line intensity). The medium 
component is centered at 99.4-99.9 kms"' (fitting uncertainty) 
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1^ 




s 



80 100 120 

Vlsr (km/s) 

Fig. 1. HIFI spectra of H'^O and H'^^O lines (in black), with 
continuum. The best model fits are shown as red lines over 
the spectra. Green dashed lines are the model fits with constant 
VfHr= 2.5 km s"' and V,„/a/;= 0.0 km s"^ . Vertical dotted lines in- 
dicate the Vlsr at 99.4 km s ' . The spectra have been smoothed 
to 0.2 km s ' , and the continuum divided by a factor of 2. 




80 100 120 

Vlsr (km/s) 

Fig. 2. HIFI spectra of Hj^O lines (in black), with continuum. 
The best model fits are shown as red lines over the spectra. Green 
dashed lines are the model fits with constant V,„r= 2.5 km s^'and 
V,„/a//= 0.0 kms"'. Vertical dotted lines indicate the Vlsr at 
99.4 kms The spectra have been smoothed to 0.2 kms ', and 
the continuum divided by a factor of 2. 
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80 



100 

Vlsr (km/s) 



120 



Fig. 3. Continuum subtracted HIFI spectra of the C'^O 9-8 and 
CS 11-10 lines (in black). The model fits are shown as red 
lines over the spectra. Vertical dotted lines indicate the Vlsr at 
99.4 km s"' .The spectra have been smoothed to 0.2 km s"' . The 
unidentified fine in absorption seen at 88.6 km s ' is indicated. 



while the broad one is redshifted by a few km s 'from the source 
velocity. Main water line profiles are more complex than for the 
rare isotopologues. Except for the 0-H2O 22i-2i2, P-H2O 2n- 
2o2, and 0-H2O 3i2-3o3 lines, none of the spectra exhibit pure 
absorption or emission line profiles. 

For low-mass objects, Kristensen et al. ( 2010| l interpret the 
naiTow emission as the signature of the passively heated enve- 
lope while the medium component might be due to shocks, but 
on spatial scales smaller than the molecular jets; the broad emis- 
sion arises in the extended molecular outflow. 

The para and ortho ground-state lines are deeply absorbed 
but their profiles exhibit a broad outflow component too. The o- 
H2O lio-l()i and P-H2O In-Ooo line profiles are fully absorbed 
as is more or less the 0-H2O 2i2-loi line while the peak depth of 
the 0-H2O 22i-2i2 line is 1.1 K (~ 65% of the continuum level, 
see Table[2]i. Thep-H20 2o2-lii line profile is also dominated by 
an absorption dip at 1.6 K (~ 11% of the continuum), but reveals 
a strong (up to IK) and broad outflow (21 kms"') component. 
This line is actually the most sensitive to the outflow. 



above 130 K reveal an infall signature characterized by an asym- 
metrical profile with a stronger blue component compared to the 
red one (see for example the P-H2O 2ii-2o2 spectra). 



3.3. Rare water isotopologues 



Among all observed HPO 



and O 



lines, the o-H'^O 1 1 



loi and p-Hl**0 2o2-lii lines are not detected (assuming a line 



upper limits are 280 and 650 mK kms 



iOl "..u F **2 

width of 6 kms ' 

respectively). The o-Hj^^O 3i2-3o3 line is tentatively detected at 
2 cr level. 

The observed profiles of lines connecting with a ground-state 
ortho or para level are dominated by absorption from the cold 
outer envelope. This is a bit surprising as one would have ex- 
pected at least the o-Hj^O 2i2-loi to be in emission because of 
the high energy level involved. All these absorptions are centered 
at 99.6 ±0.3 km s^'and can be fit by a single Gaussian with a ve- 
locity width ranging from of 4.7 to 6.4 kms"', corresponding 
to our "medium" component. No broad component is observed. 
Line parameters are given in Table |2] 



Except for the 0-H2O 221 -2i2 line with E,ip at 194 K which 
is seen in absorption, the emission lines involving energy levels 



4. Analysis 

4.1. Line Asymmetries 

Outflows, infall, and rotation can produce very specific line pro- 
files with characteristic signatures (see Fuller et al.||2005| and 
references therein). Outflow and rotation give rise to both red 
and blue asymmetric lines while the profiles of optically thick 
lines from infalling material have stronger blue-shifted emission 
than the redshifted emission. However, in some circumstances, 
outflow or rotation could also produce a blue asymmetric line 
profile along a particular line of sight to the source. 

Compared to the C"^0 (9 - 8) optically thin line (v - 
99.4 km s '), several water lines profiles show clear asymmetry. 
Focusing our analysis on emission lines, strong asymmetries are 
detected in the P-H2O 2n-2()2 and 0-H2O 3i2-3o3 lines. Both 
profiles clearly exhibit a stronger blue component, presumably 
indicating the presence of infall. Actually, strong asymmetries 
have previously been measured in the CS line profiles (involving 
lower energy levels) towards this source by |Herpin et al.| ( [2009| l, 
and the CS 1 1-10 line observed here exhibits an infall signature 
too. 

The absorption lines (including the rare isotopologues) do 
not show clear asymmetry, but as indicated by the absorption 
peak velocities (see Table [2]l, several of them are indeed red- 
shifted (relative to the source velocity of 99.4 kms"'), as ex- 
pected in case of infall. 

The P-H2O 2o2-l 11 line profile, with a very sightly stronger 
blue peak than the red one, exhibits a strong self-absorption 
dip at the source velocity. This almost symmetric double-hom 
profile might be produced by the outflow (Fuller et al. 2005 | l. 
Actually, the absorption is seen against both the outflow and 
continuum emission. Hence, as shown by Kristensen et al. (in 
press), the outflow is also embedded (i.e., the absorption layer is 
in front of both the emitting layers). 
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Table 2. Observed line emission parameters for the detected lines, ulsr is the Une peak velocity. T,^b is the intensity of the main 
peak (i.e. of the strongest emission/absorption). Tcont is the (corrected) continuum level measured from the corresponding WBS 
spectra. Ay are the velocity full width at half-maximum (FWHM) of the narrow, medium and broad components. 



Line 


"LSR 


Tmb 










T 


F 






[km/s] 


[K] 


[K] 


[km/s] 


[km/s] 


[km/s] 




[10"-' W.cm"-] 






o-Hf O lio-loi 


100.8 ±0.5 


0.40* ±0.05 


0.50±0.05 


- 


6±1 


- 


0.20 ± 0.03 








o-Hf O 3i9-3oi" 


101.2±0.2 


2.6*±0.4 


2.6±0.4 


- 


- 


- 


0.02 ± 0.01 








p-Hf O l„-Ooo 


99.8±0.3 


1.8''±0.3 


2.6±0.4 




5.8±0.3 




0.33 ± 0.07 








p-Hf O l„-Ooo 


99.6±0.1 


2.4''±0.4 


2.6 ±0.4 




6.4±0.2 




0.08 ± 0.02 








o-Hf O 2,2-loi 


98.6±0.2 


2.6*±0.5 


3.1±0.6 




4.7±0.3 




0.14 ±0.05 








0-H2O lio-loi 


99.3±0.2 


0.0* ± 0.05 


0.50±0.05 




7.5±0.2 


35.5±0.2 


> 6 


6.0 ± 0.6 


0.057 ± 


0.006 


P-H2O 2u-2o2 


97.4±0.1 


1.8±0.2 


1.0±0.1 




7.5±0.2 


18.5±0.4 




92 ±9 


0.87 ± 


0.09 


P-H2O 2o2-lll 


94.8+0.3 


1.6''±0.2 


1.8±0.3 






21.0±0.5 


0.14 ±0.01 


34 ±5 


0.32 ± 


0.05 


O-H2O 3i2-3o3 


98.0±0.2 


0.6±0.1 


2.6±0.4 




5.9±0.2 






7 ± 1 


0.064 ± 


0.009 


P-H2O lu-Ooo 


97.8±0.2 


0.0*±0.02 


2.6±0.4 




10.3±0.2 


30.5±0.6 


> 6 


8± 1 


0.07 ± 


0.01 


O-H2O 221 -2i2 


101.1±0.2 


1.1*±0.3 


3.0±0.6 




6.8±0.2 




1.0 ±0.1 








O-H2O 2i2-loi 


97.9±0.2 


0.0''±0.3 


3.0±0.6 




11.8±0.4 


20.5±0.6 


5.72 ± 0.02 


8 ± 2 


0.08 ± 


0.02 


CS 11 - 10 


98.3±0.3 


0.42±0.04 


0.50±0.05 


2.9±0.9 


10±1 












C"*09 - 8 


99.4±0.3 


0.30±0.05 


1.8±0.3 




5.6±0.3 













Notes. tentative detection. * Line in absorption. Temperature for the absorption dip with continuum 



4.2. Opacities 

For the lines in absorption, we estimate the opacities at the max- 
imum of absorption from the line-to-continuum ratio in Table [2] 
using the formula 



(1) 



2007 1 to derive the envelope temperature and density profiles. 
The physical model consists of a spherical envelope with no disk 
and no cavity to emulate a 1 -dimensional model. The density 
decreases exponentially with radius and it follows the relation: 



n(r) - no 



(2) 



and assuming that the continuum is completely covered by the 
absorbing layer. 

The Hj^^O lines are somewhat optically thick, with r a: 0.2- 
0.3, if we exclude the tentatively detected o-Hj'^O 3i2-3o3 line. 
In contrast, the H^O lines are close to optically thin, with t » 
0.1. 

Of the Hj^O lines detected in absorption, all but the P-H2O 
2o2- 1 1 1 line are highly optically thick. Note that the optical depth 
of this line may be (much) higher if not just continuum but also 
line emission is absorbed (see van der Tak et al. in preparation). 
The 0-H2O lio-loi, P-H2O lii-Ooo and 0-H2O 2i2-l()i line pro- 
files show absorption nearly down to the zero-temperature level, 
leading to opacities larger than 3 in the coolest shells. All these 
absorptions are due to cold material in which the protostar is 
embedded. 



5. Modeling 

The preceding sections have demonstrated how the line profiles 
can be decomposed into various dynamical components empiri- 
cally and how their asymmetry can signal the presence of infall. 
In this section, we model the full line profiles in a single spher- 
ically symmetric model with different kinematical components 
due to turbulence, infall, and outflow. 



5.1. Metliod 

We use the 2-dimensional Whitney-Robitaille radiative transfer 
code (hereafter WR, [Whitney et al.|2003t|Robitaille et al.|2006| 



where = 5.77 X 10^ cm"-' is the number density at a ra- 
dius ro - 1 X 10^ AU for the derived envelope mass. We as- 
sume a p value of -1.5 (model parameters are listed in Table 
3\. Dust opacities (with thin ice coated grains) are taken from 
Ossenkopf & Henning| (fT994|. The envelope temperature and 
mass are constrained by comparing the WR derived spectral en- 
ergy distribution (SED) with existing sub-millimeter continuum 
fluxes (seeFig.|5]l from our Herschel observations and |B ally et al.| 
(^0 T0|, CSO/SHA RC and IRAM/MAMBO observations from 
Mot te etal.| ( |2003] l, JCMT/SCUBA and SPI TZER-MIPS archiv e 
data (see Fig. |5J, and KAO observations ( [Lester et al.||1985| . 
Observed peak fluxes were normalized to the source size (55000 
AU ^ 10" at 5.5 kpc) derived from 1 .2 mm continuum at 3 sigma 
over noise level (Motte et al. 2003). 

Line emission is modeled with the ID-radiative transfer code 
RATRAN ( [Hogerheijde & van der""^[2000) l using the WR 
derived envelope temperature and density profiles (see Fig|4]i, 
following the method described in [Marseille et al.] (2008) and 
'Chavarria et al. (2010). The H2O coUisional rate coefficients are 
from |Faure et al., (2007 ) . 

Our model of W43-MM1 has two components: an out- 
flow and the proto-stellar envelope. The outflow parameters 
(intensity, velocity width) are derived from the Gaussian fit- 
ting of the emission lines, with FWHM between 10.2 and 35.5 
km s"' depending on the line. We have used RADEX (van der 
Tak et al.|2007| l to estimate the H2O column density and opac 



ity in the outflow. To retrieve the observed intensities for the 
outflow component we adopt Ti^,,, - 100 K, n{H2) = 10^ cm"-* 
(following the density profile shown on Fig.|4|, and a water col- 
umn density of 10'"* cm"-. Because RATRAN only enables to 



6 



F. Herpin, et al.: The massive protostar W43-MM1 as seen by Herschel-HIFI water spectra 



^ 10 




10 



10-^ 10" 
Radius [AU] 



10 



Fig. 4. Density and temperature profile derived from SED mod- 
eling. 



include the outflow as an additional component (with T, opacity, 
and width) not part of the radiative transfer process, the outflow 
component is hence assumed to be an unrelated component, both 
spatially and in velocity. The result is directly incorporated into 
the ray-tracing part of the RATRAN code (see Hogerheijde & 
van der Tak 2000). As a consequence, we stress that the absorp- 
tion in the modeled line profile is very likely overestimated, i.e.. 



the cold layer in size or temperature. 

The envelope contribution is parametrized with three input 
variables: water abundance (xhio), turbulent velocity (V,„r), and 
infall velocity (Vinf„ii)- The width of the line is adjusted by vary- 
ing Yiur- The line asymmetry is reproduced by the infall velocity 
y infall- The line intensity is best fitted by adjusting a combina- 
tion of the abundance, V,,,,., and Vouifiow parameters. We adopt 
the following standard abundance ratio s for all the lines: 450 



for Hi^O/Hi^O, 4.5 for Hi^O/Hi^O (jWilson & Rood 



1994 



Thomas & Fuller|2008| ), and 3 for ortho/para-H2U. I he rnodels 



assume a jump in the abundance in the inner envelope at 100 K 
(see Sect. |5.4| i. Table|4]gives the parameters used in the models. 

To fit the observed lines we first modeled the rare isotopo- 
logues (Hj^O, then H2'^0) lines and the Hj^O line involving the 
highest energy level (0-H2O 3i2-3()3; optically thin), using the 
same abundances for all lines (Hj^O abundances are derived 
from Hj^O and Hj^O values times the isotopic abundance ra- 
tios). Once we are able to reproduce the main features of the 
profiles by minimizing residuals in a grid of values, we model 
the rest of the lines using the same parameters, including the 
outflow component when justified. 

We simultaneously use two different criteria to quantify the 
quality of our model: 

- we quantify the error e relative to the intensity T and width 
Av of the line 



1 

"^2 



- peak,mod ^ peak,ohs 



' peak,obs 



1 

^2 



Av, 



wd 



Av 



obs 



Av, 



obs 



(3) 



Table 3. Parameters used in Whitney-Robitaille continuum 
model and derived physical quantities. 



Parameter used 


Derived parameter 


Inner radius (AU)" 


100 


Mass (Mq) 


3520 


Outer radius (AU)" 


55000 


Tou, (K) 


19 


Distance (kpc) 


5.5 


7-,,, (K) 


1377 


Luminosity (Lq ) 


2.3x 10* 






Density exponent 


-1.5 







Notes. " IMotteetal.lpOOS) 



10" 



I 



10 



10" 



10" 



-7 



10 



10 



-10 



-11 



T] 1 I 

3520 Mo 
2.3e4 
' 55000 AU 
5.5 kpc 



— I — I — I I I I I 



— 1 — I I I I 1 1 1 r- 



W43 MMl 




10 10 10 

Wavelength [/xm] 



Fig. 5. Spectral energy distribution for W43-MM1 calculated us- 
ing the Whitney-Robitaille model for the parameters in Table [3] 
EiTor bars are from the references cited in Sect. 15.11 



the similarity between the observed and the modeled line 
profiles is quantified by comparing for each channel the ob- 
served iji) and modeled intensities (T*) (above the 3cr de- 
tection limit). We define a detection function D, and a corre- 
lation function C, as follows 



Dj — { \ ^' ^ and C, 
' ' or ' 



1 if|r, -r,*| <3(r, 

or 



(4) 



The computed parameter is called the profile similarity factor 



z = 



(5) 



Values of s and S are respectively minimized and maximized by 
varying the parameters in Table |4] For consistency purpose with 
other works, we have also checked the classical values. 

5.2. Velocity structure 

In order to constrain the source dynamics, we try two different 
approaches: 1) a model with constant turbulent velocity and zero 
infall (Vtiir- 2.5 kms^'and Vm/aii- 0.0 kms"') for all lines and 
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all radii; 2) a model in which the velocity parameters vary with 
radius. Inspection of the line profiles (see Sect.|3]l shows that the 
width of the velocity components is not the same for all lines, 
nor is the infall contribution, so we do not expect a model with 
equal velocity parameters for all lines to fit the data well. Indeed, 
for most of the lines, the model 1 (see Figs. [T]and|2]) fails to 
reproduce the global shape of the profiles, especially those of 

The best model has turbulence and infall both varying with 
the radius (see Figs. [T] and |2]i, as already seen by Caselli & 
|Mye rs ( 1995 ) for the turbulence. We first tried to apply a power- 
law variation for both parameters but without any success. We 
also unsuccessfully tried a "two-steps" turbulence (Vrur- 2.2 
kms""' up to a certain radius, i.e., around 10"^ AU, and 3.5 
km s""' beyond). A "multiple steps profile" (Fig|6]l for the turbu- 
lence proved to be a much better approach to fit all line profiles 
with the same model. The infall is also best constrained by a step 
function given in Fig|6] Each shell of the envelope (at radius R) is 
hence characterized by a temperature, a density, a turbulent, and 
an infall velocity as shown in Fig|4]and|6] We have estimated the 
uncertainties of these velocities by making the velocity varying 
around the "best value " at each radius. 

Line profiles for all species are well reproduced by our model 
(Figs.[T]|2]and|3]l, with E similarity coefficient larger than 90%, 
and minimized s andx^ for most lines (o-Hj'^O lio-loi, P-H2O 
2n-2o2, and 0-H2O lio-loi have E ^ 75%), even with an out- 
flow contribution fit with a simple Gaussian model. The model 
predictions are consistent with the upper limits for the non- 
detected lines o-H^O lio-loi and p-Hj^O 2()2-lii, and the ten- 
tatively detected o-H^^O 3i2-3o3 line. The CO and CS lines are 
well fit by the model derived from the other lines. All modeled 
lines are centered on Wlsr equal to 99.4 kms but the outflow 
component central velocity varies somewhat from line to line 
by ±0.5 - 1.0 kms"'. In addition, we apply our model to the 
Hi'^O 3i3 - 220 line at 203.3916 GHz (£„p=204 K) observed 
with t he IRA M-30m telescope (AV = 3.8 kms"') by Marseille 
et al. ( 2010a| l. The fine profile is perfectly reproduced. 



Table 4. Physical quantities derived from the RATRAN model. 



5.3. Infall and accretion rate 



Actually, outflow (e.g.. Walker et al.||1988 1 and rotation ( Zhou 



1995 Fuller et aL]|2005[ l can also give rise to blue asymmetric 



lines along a particular line of sight to the source. But consider 
ing that blue asymmetry is observed in several water transitions 
(see also Herpin et al. 2009 for CS), and that there is no evidence 
of contradictory information in the form of red asymmetric pro- 
files, the infall explanation is the most likely. Nevertheless, only 
mapping (and position-velocity diagram study) with a high an- 
gular resolution (of a few 0. 1 arcsec) could help to separate the 
effects of rotation and infall. Indeed, rotation might produce a 
symmetry reversal in outer part in direction perpendicular to the 
rotational axis where the optical thickness is less ( Zhou|[T995 
|Yamadaetal.|2009| ). 

From the infall velocity Vi„faii we estimate the mass accre- 
tion rate Mace applying the method described by 'Beltran et al. 
( 2006| l. As a first rough estimate we assume that the gas is un- 
dergoing spherically symmetric free-fall collapse onto the cen- 
tral object. For any radius R, the mass accretion rate is given by; 



Mr,rr - 47TR^mn Vi 



infalU 



(6) 



where m is the mean molecular mass, n the gas volume density, 
and Vinfaii the infall velocity at R as shown on Fig. |6] 



Parameter 



%20 

Post-jump XH20 
o/p 

^180/"0 

VturCkms-') 
^outflow (km s"') 
Vinfall,max (kms" 
Vlsr (kms-') 



8.0 (±1.0)xlO-*' 
1.4 (±0.4) X 10"" 
3 ±0.2 
4.5 
450 
2.2-3.5 
10.2-35.5 
-2.9 
99.4 



The maximum infall velocity of 2.9 kms"' is mainly de- 
rived from the P-H2O 2ii-2o2 and 0-H2O 3i2-3o3 line profiles. 
From the model (Fig. |6] and Sect. 5.2 1 we can estimate for that 
infall velocity the corresponding radius from the protostar to 
5.7-7.3 X 10'^ cm (38100-48800 AU). The inferred mass accre- 
tion rate for this range of radii is then 3.5 to 4.0 xlO"^ Moyr~' 



slightly lower by a factor 2 than what is derived by [Herpin 
|et al.| ( |2009| ) from CS observations. At such radii, we are ac- 
tually likely measuring the accretion from the envelope toward 
the disk, if there is a disk. However, according to our model, the 
infall velocity is varying with radius, i.e., the mass accretion rate 
is not constant with radius. Inside of 26800 AU, the infall ve- 
locity is only 0.4 kms"', which leads to a mass accretion rate 
between 4.0 xlO"" and 5.0 xlO"^ Moyr"'. 



5.4. Abundance structure 

Modeling the whole set of observed lines allows us to determine 
the abundance structure. Among the lines mostly excited in the 
inner envelope, only the higher-7 Hj^O (o-Hj^O 3i2-3()3 andp- 
ii\^0 3i3-32o) and the 0-H2O 3i2-3o3 are optically thin enough 
to probe the inner part of the envelope. An abundance jump at 
100 K is necessary to reproduce these lines while the other lines 
are not sensitive to it. The derived Hj^O abundance relative to 
H2 is 1.4 xlO""* in the warm part where T >100 K while it is 
8.0 xlO"^ in the out er part where T <100 K , in agreement with 



Marseille et al. 



(|2010all from Hf O obser- 



the values found by 

vations. No deviation from the standard o/p ratio of 3 at high 
temperature (>100 K) is found. Hence, the one-dimensional ap- 
proach seems acceptable as the shape of all line profiles is well 
recovered. 

From the integrated fluxes given in Table|2](measured for the 
lines at least partly in emission), one can derive the water lumi- 
nosities. We then estimate, assuming isotropic radiation and at a 
distance of 5.5 kpc, that the minimum total HIFI water luminos- 
ity is 1.5 L0 (sum of all individual observed luminosities). The 
true water emission from the inner part could be much greater 
but the cool envelope absorbs much of the emission. Moreover, 
from the modeling we estimate the total water (radial) column 
density in the envelope to be 1.56 10^' cm"^, corresponding to 
0.11 Mo (M„i8o = 2.4 X 10"^ Mo and M„no = 5.6 x 10"^ Mo). 
More than 97% of this mass is in the inner part. The CS and 
C'**0 column densities are respectively 7.3 x 10'^ and 8.2 x 10'^ 
cm"^ leading to masses of 3.2 x 10""* (CS) and l.Ox 10"^ (C"^0, 
5.0 for CO) Mo. Hence the mass of oxygen trapped in H2O and 
CO is around 3 Mo, 96% locked in CO and 4% in H2O. 
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5.5. Uncertainties and robustness of ttie parameters 

Uncertainties in our calculations result mainly from three 
sources: 

- the intensity calibration of the observed lines (see Sect.[2|. 

- The physical limits of our model whose radiative transfer 
computations are based on a physical description of the 
source derived from a ID SED modeling. 

- The parameter sensitivity of the line radiative transfer model. 

We have reiTin the models first increasing, then decreasing 
the observed line intensities by the uncertainties given in Section 
121 We found differences in the estimated abundances between 
10 and 30% depending on the line. To test the influence of the 
physical model, we decreased the envelope mass by 30% (which 
gives an SED just below the error limits) and recalculate the wa- 
ter abundances. We found diff'erences up to 75%. In addition, 
considering that the error on the exponent of the density pro- 
file is 0. 1 (see Sect. |5.1| l, we have tested the robustness of our 
results when varying the density profile within that range: we 
found that the abundances are impacted by 30-40%. Moreover, 
the uncertainty for the distance to the source (we adopt 5.5 kpc. 



E„p=248 K 



137 114 61 53 



but Nguyen Luong et al. 



2011 



place the distance at 5.9;^^^ kpc) 
translates mto an uncertamty on the physical model: 10% ifor the 
distance leads to 15 % for the mass hence 30 % in the abundance. 
This of course impacts our results. In addition, we checked that 
varying the best fit parameters from Table |4]by 5 % has no influ- 
ence on the derived abundances. Variation of the turbulent veloc- 
ity within the error bars (and varying the radius by 15%) shown 
on Fig|6]does not affect the results. Actually this turbulent ve- 
locity, as explained in Sect. |5.2[ is well constrained thanks to the 
high number of observed lines, efficiently probing different ra- 
dius. The determination of the infall velocity versus the radius is 
more uncertain because the infall signature is not visible for the 
quite completely absorbed lines. 

Unfortunately, we are not able to determine if these un- 
certainties are independent of each other. Hence, only the un- 
certainties coming from the HIFI intensity calibration provide 
us with a realistic quantitative estimate of the impact on our 
calculations. In a forthcoming study, we plan to investigate 
how a 2-dimensional model incorporating disk and cavity (e.g., 
[Hosokawa et al.|2010| l would impact the results presented here. 
The fact that the spherical model provides such a good fit argues 
against large changes. Nevertheless the modeling results pre- 
sented here are the best we got, but it is obvious that the model 
is very likely not the only possible one. However, we stress that 
by varying parameters, the derived numbers change, but not our 
main conclusion that the turbulence and the accretion rate vary 
with radius. 




R (AU) 

Fig. 6. Variation of the turbulent and infall velocities, respec- 
tively Viurh and Vi„f, with the radius R (AU) to the central object. 
Overplotted are the radii corresponding to the temperature jump 
(dashed) in the model and to the radius where lines of different 
energies are mainly coming from with upper energy levels as 
dashed lines. 



the rare isotopologues mainly comes from the medium compo- 
nent in contrast with the low-mass objects where the broad out- 
flow is the main origin. Second, no narrow component is seen 
in the water lines, maybe because of opacity effects. The nar- 
row component emission comes very likely from the hot core 
neighborhood (i.e., the passively heated inner envelope) and ob- 
servations of higher-energy level water lines are needed to probe 
more deeply this region. The broad component originates in 
the molecular outflow similar to what is observed in low-mass 
objects while the medium velocity component is likely due to 
combination of turbulence and infall. The main difference with 
W3IRS5 is the presence of infall and no sign of expansion. 

According to our model, the outer water abundance derived 
in W43-MM1 (8 x 10"**) is several times higher than what is 
found in the mid-IR bri ght HMPO W3IRS5 (j^ou, = 1.8 x 10 ^ 
Chavarria et aL]|2010| or in the HMC G31.41 (^out = 10 ** 
Marseille et al.||2010b i. It is also an or der of magnitude large r 
than in low-mass objects (xou, = 10"^, Kristensen et al.|2010 l 



In contrast, although observing lines with higher upper energy 
levels could help to improve the estimate, the inner abundance 
in W43-MM1 (1.4 x 10"^) can be constrained to a value sim- 
ilar to what is f ound in other h igh-mass sources (e.g. = 



1.0- 2.0 X 10-4 
i.e. , most oxygen contained in water. 



van der Tak et al. 2006; Chavarria et al. 2010 1, 



6. Discussion 

6. 1 . Comparison witli previous work 

These observations towards W43-MM1 of a large set of water 
lines, together with CS and C'**0, reveal that, as for low- and 
intermediate-mass obj ects (^K ristensen et al. 2010; Yild iz et al.| 
[Johnstone et al.pOlO) , one finds a broad outflow (> 20 



|20T0l p 



kms"') and a medium (5-10 kms ') velocity component. On 
the opposite, the naiTow (< 5 kms ') component is not clearly 
detected, even if it might be present in the CS line profile. They 
show two features not seen in lower mass types of objects. First, 
as noted by Chavarria et al. (2010 1 in W3IRS5, the emission of 



6.2. W43-MM1 and tine ir\igir\-mass star formation process 

This work shows that the observed molecular emission is dom- 
inated by turbulent velocities (v,„r > 2.2 kms"') greater than 
the speed of sound (as - 0.3-2 kms"' at 20-1000 K) at 
all radii. This may be the result from the fact that our model 
is 1 -dimensional, or from the possible role of rotation in the 
line broadening. We are not able to test the latter effect with- 
out mapping W43-MM1 at higher spatial resolution. However 
large rotational velocities (> 10 km s"') are likely to be found at 
radii closer than 600 AU ( Keto & Zhang|[20T0) , and moreover 
its main impact would be a small velocity shift of the molecular 
line profiles ( iBelloche et al.|2002| . This high level of turbulence. 
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derived from so many indicators, agrees with the conclusion of 
Plume et al. ( 19971) that regimes of massive-star formation are 



highly turbulent. Figure |6] illustrates how the turbulent veloc- 
ity varies with distance to the center of the object as derived 
from our model. We can then state that the turbulent motions in- 
crease with radius, in agreement with the turbulent core model 
of|McKee & Tan| (12003 ) and Larson's studies in clouds ( [Larsonl 
2003|. The lower degree of turbulence in the inner envelope can 



be due to the rotation and the density structure suppressing tur- 
bulence close to the protostar 

Additionally, Fig. |6] shows the radii where the lines are 
mostly excited in our model (lines are dominated by one main 
turbulent velocity): lines with low Eup are more sensitive to outer 
regions (i.e., higher energy level transitions come from the inner 
part). 

Moreover, the free-fall accretion rate would be around 6.4 x 
10"^ Moyr"' (for as - 0.3km s"') while here we detect ultra- 
sonic gas motions and mass accretion rate between lO"** and 
10"^ Mgyr"' depending on the radius. Nevertheless, we stress 
that this very high derived accretion rate assumes spherical ac- 
cretion, which may not be true. Also, if there is more than one 
central object in W43-MM1, the accretion rate per star could 
be considerably lower We can estimate the corresponding ac- 
cretion luminosity, for the largest Mace we derived, applying the 
following formula: 



(7) 



Considering a mass of 20 Mq within a radius of 100 Rq for the 
protostar ( Hosokawa et al. 20 1 ), we get an accretion luminosity 
of the order of 3 X 10"* Lp . Comp aring this value with the stellar 



luminosity of 10^ Lq derived by [Hosokawa et al.| ( |2010 ) for a 
protostar with an accretion rate of only 10"^ Moyr"^ and 20 Mq, 
we conclude that the derived accretion luminosity is unrealistic, 
compared to the observed total luminosity. As a consequence, 
the assumption o f a sing le object might not be correct. However, 
[Hosokawa et al.| pOlO) ) also show that realistic accretion rates 
are much lower (by about one order of magnitude) than the val- 
ues obtained from the simple formula we use, and as a conse- 
quence M might be more likely of the order of 10"^ Mgyr"' and 
the accretion luminosity of the order of a few 10^ Lq , hence con- 
sistent with the observed total (stellar -naccretion) luminosity. 

The derived accretion rate, although uncertain, is high 
enough (larger than 10 Moyr ') to overcome the radiation 
pressure due to the star luminosity ([McKee & Tan|2003 Yorke 
[& Bodenheimer[|2008) [Hosokawa & Omukai[|2009| ). However, 
the protostar cannot reach the zero age main sequence by steady 
mass accretion if the accretion rate is so high ( [Hosokawa et al. 
[2010| ). Of course, one does not know if W43-MM1 wiU maintain 
such a high accretion rate over a long period of time. 

As recalled in the introduction, W43-MM1, though a very 
young source from indicators such as the SED, has already de- 
veloped a hot core which is very likely at the origin of the high 



abundances of S-bearing species like H2S and OCS ( Herpin et al 



2009| ) whose evaporation from grain surfaces could 



->e important 

at such temperatures. Also, Herpin et al. (j2009) underlined the 
possibility of strong shocks occurring in the inner part (where 
T >100 K). Moreover, our study shows that the water abun- 
dance in W43-MM1 (in the outer envelope) is larger than in more 
evolved massive objects such as W31RS5 and G31.41, which 
suggests fast desorption from grains. Hence, W43-MM1 seems 
not to fit in the assumed evolutionary diagram. Actually, there 
is some indication that high water abundances are not clearly 



linked to luminosity, mass, temperature or assumed evolution- 
ary stage of the sources ( [Marseille et al. 2010b ). On the con- 
trary, the high water abundance might be related to the presence 
of a HMC in the object. But, if so, one has to explain how high 
temperatures in the inner region can boost water production in 
the outer region. Another option for high water production in 
the inner regions is by shocks. As revealed by the violent gas 
motions detected in the envelope, W43-MM1 exhibits a huge 
infall and large turbulent velocities. Infall is generally associ- 
ated with strong and fast outflows (e.g. Klaassen et al.|20lT Liu 
jet al. 201 1\ . As a consequence, rather than being due to the pres- 
ence of a HMC, the higher water abundance might be due to the 
presence of fast gas motions (i.e., shocks) as seen in low-mass 
objec ts ( [Codella et al.|2010[[Nisini et al.|20T0[|Kristensen et al. 
'2OIOI) where H2''0 abundances up to 10 ^ can be found along 
outflow walls. 



7. Conclusions 

We have presented here Herschel-HIFI observations of W43- 
MMl, for the first time, of as many as fourteen far-IR water lines 
(H'^'O, H'^O, Hi**0), CS(ll-lO), and 01*^0(9-8) fines. We first 
analyzed the line profiles, then we modeled the observations us- 
ing a radiative transfer code. We estimated outflow and infall 
velocities, turbulent velocity, and molecular abundances. 

Compared to what is observed towards low- and 
intermediate-mass objects, the water lines observed towards 
W43-MM1 are obviously brighter, but reveal the same broad 
and medium velocity components. The narrow component, 
from the hot core seen in W43-MM1, is not so commonly 
detected, likely because of opacity reasons. Hence observations 
of higher-energy water transitions are needed to probe this hot 
core. 

From the modeling we estimate an outer water abundance of 
8 10"*^. This value is higher than observed in other sources and 
could be due to the fast outflow and infall observed in W43- 
MMl. Moreover, the high level of turbulence derived in this 
source seems to be the usual regime in massive protostars, while 
the huge infall, i.e., mass accretion rate, very likely makes W43- 
MMl a special case. Both turbulence and accretion rate vary 
with radius. If confirmed, the measured high level of turbulence 
combined with a high accretion rate agrees with the turbulent 
core model, although we cannot yet rule out without high spa- 
tial resolution data the competitive accretion model. Finally, the 
estimated accretion luminosity is high enough to overcome the 
expected radiation pressure. 

HIFI mapping of the outflow, to be presented in a next paper, 
will help to determine the detailed kinematics of this compo- 
nent. In addition, the gas cooling budget will be studied in detail 
through a complete census of aU water, plus [OI] and OH fines 
with PACS. 
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Appendix A: Absorption lines on the line-of-sight: 
cold foreground gas ? 

The ground-state ortho- and para-Hj^O and the 0-H2O 2i2- 



Ifli spectra (see Fig A. 1 1 reveal the presence of deep absorptions 
at several velocities (67.2, 71.0, 79.5, 82.5, and 87.8 kms"'; de- 
rived from Gaussian fits). The parameters of the absorption fea- 
tures are listed in Table A. 1 For each velocity component we 



derive the optical depth at the maximum absorption dip, based 
on the assumption that the excitation temperature is negligible 
with respect to the continuum brightness temperature. The op- 
tically thickest absorption is at 79.5 kms"' (ti669 ^ 3.1), fol- 
lowed by components at 67.2 and 82.5 kms"' while the two 
other absorptions are optically thinner. Because the absorption 
of the continuum signal at 79.5 kms ' is almost complete, the 
absorbing material must be cold. The Gaussian fits reveal that the 
67.2, 71.0 and 87.8 kms"' absorptions are 1.4-1.9 kms"' wide 
while the absorption at 79.5 km s"' is definitely broader (5.3-5.9 
kms"'). 

Assuming that all water molecules are in the ortho- and para 
ground states and that the excitation temperature is negligible 
with respect to the continuum brightness temperature, we cal- 
culate the column densities from the opacities applying the 
following formula: 



(A.l) 



where t is the optical depth, v the line frequency, Sv the 
width, c the speed of light, A the Einstein coefficient, and gi and 
gu the degeneracy of the lower and upper levels of the transition, 
respectively. 

Low column densities are derived, of the order of a few 10" 
cm"^, except for the 79.5 km s"' component where is several 
times larger (these values of the column density explain why the 
absorption are not seen for other transitions). From these calcu- 
lations we derive the o/p ratio for each component. Surprisingly, 
if the absorptions are due to cold foreground clouds, the o/p ra- 
tio does not strongly diverge from 3 if we take into account the 
uncertainties. This means that, in these clouds, water is produced 
in what is called the high temperature regime (Tjp,„ > 50K). Of 
course, we stress that opacities might be underestimated because 
of the assumptions we made. 
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Table A.l. Opacity, total column density and o/p ratio derived for each line-of-sight absorption component detected. 
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Fig. A.l. Overplotted spectra (divided by the continuum) of the 
ground-state ortho- (black) and para-Hj*"© (red) and the 0-H2O 
2i2-loi(blue) lines. The different absorption components, apart 
from the main Hlf'O emission from the object, are shown by a 
dashed line and the corresponding velocity is given. The spectra 
have been smoothed to 0.3 kms"'. 
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